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Abstract 

We reconsider the applicability of classical nucleation theory (CNT) to the calculation of the 
free energy of solid cluster formation in a liquid and its use to the evaluation of interface free 
energies from nucleation barriers. Using two different freezing transitions (hard spheres and NaCl) 
as test cases, we first observe that the interface-free-energy estimates based on CNT are generally 
in error. As successive refinements of nucleation-barrier theory, we consider corrections due to a 
non-sharp solid-liquid interface and to a non-spherical cluster shape. Extensive calculations for the 
Ising model show that corrections due to a non-sharp and thermally fluctuating interface account 
for the barrier shape with excellent accuracy. The experimental solid nucleation rates that are 
measured in colloids are better accounted for by these non-CNT terms, whose effect appears to be 
crucial in the interpretation of data and in the extraction of the interface tension from them. 

PACS numbers: 64.60.qe, 68.03.Cd, 68.35.Md 
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The decay of metastable states, such as the solidification of a supercooled liquid, takes 
place through the nucleation and growth of some small-sized droplet within the system 4]. 
The initial stage of the phase transformation is usually described within the time-honored 
classical nucleation theory (CNT) js — 7 1 , where the droplet is envisaged as a sphere of, say, 
bulk solid, separated from the liquid by a sharp interface, giving rise to a free-energy penalty 
proportional to the interface area and a total Gibbs-free-energy activation barrier 

AG(n) = -\Afi\n + An 2/3 , (1) 

where n is the number of particles in the solid cluster, A/i < is the chemical potential 
difference between solid and liquid, A = (367i~) 1//3 p s 2//3 cr with p s the bulk-solid number density 
and a the specific surface energy (surface tension) of the planar interface, all anisotropies 
being neglected at this stage. The droplet grows if it exceeds a critical size n* corresponding 
to the maximum AG(n) (= AG*). CNT is routinely used to estimate the nucleation rate 
/ = Joe _/3AG ", where (3 = 1/(&bT) is the inverse temperature and 7n a kinetic prefactor that 
varies slowly with T. Clearly, this connection between / and a [7] relies on several severe 
approximations. First of all the choice of an appropriate reaction coordinate, here the 



droplet size n, an issue largely discussed and criticized in the literature [8l-lllj. Moreover, 
Iq is notoriously influenced by genuinely non-equilibrium effects and various expressions 
resulting from a more detailed consideration of the nucleation kinetics are known since a 



long time 



3j4i4|. 



In this Letter, we do not address the issue of the validity of CNT for predicting the 
nucleation rate but rather consider an even more fundamental question, namely the efficacy 
of CNT in describing the dependence of the interface free energy of the solid cluster on 
its size. Our starting point is to show that the profiles of AG(n) obtained by numerical 
simulation of nucleation clusters in a variety of systems are not consistent with Eq. (CD). We 
then explore corrections, some already present in the literature, some novel. It emerges that 
the numerical profiles can be accurately reproduced by assuming a diffuse and thermally 
fluctuating solid-liquid interface. Finally, we show how this finding is of direct use to inter- 
pret nucleation rates and correctly extract interface free energies from them, a result that 
should be of considerable interest to experimentalists. 



We begin by displaying in Fig. 1 existing accurate simulation data 



cluster nucleating inside a bulk liquid, available for hard spheres 



or AG(n) of a solid 
151 ] and for the Fumi- 
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FIG. 1: (Color online). Gibbs-free-energy cost AG(n) of an n-cluster in units of kgT. Left: hard 
spheres at a packing fraction of 0.5207, from [15]; right: NaCl at T = 825 K, from jig ]. Top panels: 
blue crosses, MC data (selection of data points); red dotted line, CNT best fit; black solid line, 
Landau-theory best fit. In all fits, data for n < n*/5 are ignored. Bottom panels: deviation of the 
fitting curves from the data. CNT (S = e = C = in Eq. ©); CNT + offset (6 = C = 0); CNT + 
Tolman (e = C = 0); Landau (C = 0). Values for a from CNT are: hard spheres, f5ad 2 = 0.724 
(sphere diameter d); NaCl (cubic nucleus), a = 79.75 erg/cm 2 . In hard spheres, the optimal f3ad 2 
is 0.737 for Landau theory (with 6 = -0.017 d and ? = -0.332 d 2 ), 0.741 for CNT + offset (with 
4irf3ae = -2.599), and 0.761 for CNT + Tolman (with S = 0.086 d). In NaCl (cubic nucleus), 
the optimal a is 88.87 erg/cm 2 for Landau theory (with 5 = 0.716 A and ? = 0.471 A 2 ), 83.74 
erg/cm 2 for CNT + offset (with 6/3ae = -3.833), and 88.53 erg/cm 2 for CNT + Tolman (with 
5 = 0.675 A). Note the large improvement over CNT obtained with just one more fitting parameter 
and how similar is the quality of the fit based on CNT + offset to the Landau fit. 

Tosi model of NaCl [if]] . In each case, we superpose a CNT least-square fit to Eq. ([T|) 
for comparison. It is clear that CNT is not generally adequate to describe AG(n). The 
deviations are systematic and of different sign at low and large n. The fit quality does not 
improve by restricting data to large clusters only, indicating that even in the barrier region 
the cluster free energy does not obey Eq. (JT]). To shed light on this failure of CNT, we relax 
the approximations leading to Eq. (JT|) one at a time. Eventually, we shall get a more general 
expression for the free-energy cost of a n-particle cluster for large n, which turns out to have 
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the Dillmann-Meier [13 form 

AG = AnR 2 a f 1 - % + J - \^Ps\^\ + Cln - , (2) 
\ it R 2 I 3 a 

with i? = pTry^Trps)] 1 / 3 and a a microscopic length, and where cr,5,e, and C are theory- 
dependent parameters. The assumption in (|2j) is that of a spherical cluster shape — a 
different shape, that would be determined by free-energy anisotropy, would slightly change 
the value of a but not the physical discussion that follows. The first modification to CNT 



we consider is dropping the sharp- interface approximation 18]. Within Landau theory, the 
free-energy cost of the critical droplet is the unstable stationary point of a functional, e.g. 

G[<j>] = J d 3 x {|(V0) 2 + ^(V 2 0) 2 + s(0(x))} . (3) 

Here c, k > 0, 0(x) is the "crystallinity" order parameter (OP) that distinguishes the solid 
(0 > 0) from the liquid (0 = 0), and g(<p) is the Landau free energy per unit volume of 
the homogeneous system. Below melting, T < T m , g shows, besides the liquid minimum 
g(0) = 0, a second and deeper solid minimum. Right at T m , we assume g(<f>) = c 2 o0 2 (l — 
0/0 s o) 2 (1 + t0/0 s o) with C20 > and r > — 1 [ji3], where <p s o is the value of in the bulk solid 
at coexistence, and where a non-zero value of r creates an asymmetry between the liquid 
and the solid minimum. We also assume that, slightly below T m and at fixed pressure, g 
acquires a linear dependence on AT = T — T m only through its 2 term, which becomes C20 2 
with c 2 = c 20 + c' 20 AT. With this standard setup, the free energy of a cluster of radius R is 
G[0 b], where 4>r(t) is the spherically-symmetric OP profile of the cluster. Assuming, as in 



20] . that for small supersaturation and large R, 4>n{r) may be approximated with <fio(r — R), 
where <J)q(z) is the OP profile for a planar interface centered at z — 0, the cluster free energy 
takes precisely the form (j2J), with C — 0, p s Au = c 2Q 2 o AT, and a, 5,7 all linear functions 
of AT, expressed in terms of c, k, and <po{z) [211 ] . At coexistence and to first order in the 
deviations from the 4 theory (viz. k = t = 0), we have for cr(T m ) = cr, etc. 



a = ^ f 1 + lr + Z^\ , 5=^t, and 
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where £ = a/2c/c 2 q is a measure of the interface width. 



4 



A second effect that is absent in CNT but present in nature and observed in simulations 
is that shapes of clusters, far from being static, fluctuate widely away from their mean 



shape 22|, |23|. To describe shape fluctuations, we employ a field theory for the Canham- 



Helfrich (CH) Hamiltonian, containing spontaneous-curvature and bending-energy terms in 
addition to surface tension. A CH interface Hamiltonian H s can be derived from the free- 
energy functional ([3]) for small deviations of the interface from planarity. Denoting by E the 



generic closed-surface profile and by n its outward normal, we obtain 



2l| 



U s = J dS ^x-a5V-n+^A(V-n) 2 ^ 



(5) 



where a and 5 are the same as in Landau theory and A = K<b\ n l (31) under the same hy- 
potheses for which Eq. (j3J) holds. AG(R) can be evaluated explicitly [21] for a quasispherical 
cluster 24j, where only quadratic deviations from sphericity are kept. The wavelength of 

— 1/3 

surface undulations is cut off at a lower limit a = p s to account for the granularity of 
matter. We find that the surface free energy has a form consistent with Eq. O), with new 
T-dependent parameters a, 5, and e (whose explicit expressions are given in 2l|]) and with 
C = — (7/3)&bT, which shows that small deviations around a nominally spherical cluster 
shape simply add a universal logarithmic correction to the mean-field functional form of AG. 
This correction is responsible for the well known R* 7 ^ 3 term in the exponential prefactor of 
the nucleation rate 251 ] . 

Clearly, the parameters in Eq. fl2]) are determined by the values of c and k in fl3]), as 
well as by the form of g(4>) — all system-dependent quantities that require a case-specific 
theory. We here aim at elucidating the relative importance of the different terms 5, <T, C 
implied by interface thickness and shape fluctuations. To get a quantitative measure of 
that, we directly fit the parameter values in (J2]) to the numerical results for AG(n) for 
the two systems of Fig. 1. Consistently with the assumptions underlying our mesoscopic 
description, each fit is made only to data points for sufficiently large n. We first include the 
leading ex dn 1 / 3 ("Tolman" 30j) correction to CNT. As shown in Fig. 1, this improves the 
quality of the fit significantly. The error is reduced substantially in both systems, although 
not monotonically. Only a marginal improvement is obtained if both 5 and (T are allowed in 
the fit. The inclusion of the logarithmic shape correction gives no further appreciable gain. 
Next, we attempted fitting the data by retaining just the offset (= inae) in (J2]) beyond 
a. Alone, the simple offset gave an improvement of about the same quality as with all 



5 



terms allowed. We conclude that corrections to CNT exclusively deriving from a fluctuating 
cluster shape appear to be much smaller than those arising, already in Landau theory, from 
allowing a non-zero thickness of the interface (shape fluctuations are not an ywa y immaterial 



since they renormalize, even significantly, the Landau-theory parameters [21]). Moreover, 
either the Tolman correction or, alternatively, the constant offset each lead to significant fit 
improvement over CNT. The origin of both terms is in the finite thickness of the interface, 
which makes the reversible work to create a cluster systematically smaller than what would 
be needed for the same cluster with a sharp spherical interface. 

The existing simulation data do not permit to assess the relative importance of the two 
smooth-interface contributions 5 and ? and of the logarithmic correction; more specific work 
is needed in order to decide that case by case. Using the 3D Ising model as a test system, 
we carried out extensive simulations at moderate supersaturations, computing the cluster 
free ener gy f or the nucleation process of magnetization reversal by the same method as in 



Refs. 



26 



29|. We computed AG(n) for a number of values of the field h (0.35, 0.40, . . . , 0.65, 
in J units) and plotted the ratio cr(n) of the surface free energy F s = AG(n) + |A//|n to 
the area (367r) 1//3 (na 3 ) 2//3 of the cluster surface as a function of the inverse radius n -1 / 3 (see 
Fig. 2). We verified that, for all h values considered, clusters close to critical indeed contain 
the vast majority of up spins in the system, coherently with the physical picture at the 
basis of our theories. It is evident that only the joint consideration of 5 and «f is able to 
reproduce the upward concavity of cr(n) as a function of n~ 1//3 in the n region (n > 40) 
where F s (n) oc n 2 ^ 3 . A positive offset e > is confirmed, as expected from Landau-theory 
results for e and from the formula for A in Eq. ([5]) 21] . In our regime of h, the logarithmic 
term does not change the quality of the fit; as shown in Fig. 2, this correction becomes 
sizeable only at values of n outside the fit range. However, inclusion of the logarithm has 
consequences on the optimal 5 values, which reduce from ps 0.10 to ~ 0.02 throughout the h 
range considered (the Tolman length 8 3o| is zero for the Ising model at coexistence 20]). 
From this example we conclude that a) neither 5 nor e, both arising from the finite interface 
width, can generally be neglected in the description of the nucleation free-energy barrier; b) 
shape fluctuations improve the description especially for small cluster sizes (n < 80). 

Far from being academic, the existence of these corrections to CNT has a direct impact 
on the understanding of experiments, in particular on the all-important extraction of the 
interface free energy a from measured nucleation rates. Assuming the standard activated 
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FIG. 2: (Color online). The cluster interface free energy of the 3D Ising model in units of J /a 2 (a 
being the cubic-lattice spacing), plotted as a function of ra -1 / 3 for various h values. The temperature 
is T = 0.6 T c ; starting at T with all spins down, the system is quenched to h. Two spins are part of 
the same cluster if there is an uninterrupted chain of up spins between them. The lattice consisted 
of 20 3 sites; for h = 0.35 a calculation on 25 3 sites led to practically the same AG(n). Three 
different instances of umbrella potential were considered and all led to the same AG(n) to within 



less than 0.1 ksT 



31]. Thick blue lines, MC data for AG(n); black lines, least-square fits of the 



n > 80 data points for h = 0.35 and 0.55, based on various extensions of the CNT (see legend). 
Note that only the full Landau expansion captures the upward curvature of AG(n), and that 
especially the shape fluctuations capture that of smaller clusters. In the inset, different ways of 
extracting the interface tension a out of finite-/i values of cr(oo) are compared: triangles, CNT; 
crosses, Landau theory (C = in Eq. fl2J)); open dots, quasispheres (C = — (7/3)&bT). Linear 
extrapolation of data points at h = yields a ~ 1.60. The black dot is the value of a calculated 



for the (001) interface (from Ref. 32]). The red shading indicates a estimates from Eq. ([6]) for 
cluster shapes intermediate between cubic (lower end) and spherical (upper end). Gratifyingly, 
the a values extrapolate as they should to a o"o which is higher than that of the (001) interface, 
and intermediate between cubic and spherical shape - the average shape being also intermediate 
between the two. 
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FIG. 3: Nucleation-rate data for the solidification of a colloidal fluid (Fig. 7b of 35f] ) , showing 
Y = ln(J/Jo) as a function of a quantity akin to X = (T m /AT) 2 (<£ crystal is the solid volume 
fraction). The dotted line gives the slope from which the surface tension should be extracted. 



expression 
(T m /AT) 2 



or /, a can be extracted from the slope of Y = ln(/// ) as a function of X 



33]. If CNT were exact, this slope would be a constant throughout the region of 



liquid metastability. When the more general Eq. ([2]) is employed for AG*, the slope depends 
on the distance from coexistence AT, as demanded by non-zero values of a, 5, e, C and their 
rates of variation with AT. Close to coexistence, one can write Y(X) as a power series in 
AT: 



ln(/// 



aT 2 , a'T„ 



+ 0(1) 



(6) 



AT 2 | AT, 

Here a takes the same value as in CNT, a = 167ro" 3 /(3A;BT m p 2 L 2 ra ) with L m the latent heat 

of melting per particle. However, a' is not universal: 

3o-'T m 3p s L m S 
a H 



n = . | I -! : I , (7) 

a a J 

taking a = a + a'\ AT| + . . . close to coexistence. For instance, for the 4 theory it turns out 
that a'/a = 1 + 3p s £L m /(2a) > 0. 

Due to ([()]), Y{X) develops a concavity, which is upward if a 1 > 0, as is the case for 
example in colloids (see below). The very important practical consequence is that the solid- 
liquid interface free energy at coexistence (a), the key quantity which one wishes to extract 
from nucleation rates, is determined by the slope (a) ofY(X) at asymptotically large X and 
not from the slope, generally different, at small X. As an example, deviations from linearity 



S 



in the Y(X) plot are experimentally evident in colloids, see e.g. Refs. [34], |35|. Fig. 3 
shows how data should be read to extract a. Since this procedure is not to our knowledge 
universally followed, this suggests that at least some tabulated a values may need a revision. 
Because the finite-interface corrections reduce the barrier height compared to CNT, it is to 
be expected that the true interface free energies are substantially smaller than believed so 
far. 
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Supplementary Material 

Santi Prestipino, Alessandro Laio, and Erio Tosatti 



I. LANDAU THEORY OF NUCLEATION 

In the following, we shall refer to the main text of the manuscript as MT. 

A Landau theory of nucleation based on a Cahn-Hilliard-like functional [l| describing the 
free-energy cost of a diffuse interface between two phases — "solid" and "liquid" , neglecting 
anisotropies — can be formulated as follows. Assume for simplicity a scalar order-parameter 
(OP) field 0(x). The Landau free energy is (Eq. (3) of MT): 

m = J ^ {|w) 2 +f(v 2 0) 2 +^(0(x))} , (i) 

where c, k > are "stiffness" parameters and g(0) is the specific free energy (i.e., Gibbs 
free energy per unit volume) of the homogeneous system, the bulk liquid being the reference 
state where = 0. Exactly at coexistence, the OP values are 0_ = 4> s q in the bulk solid and 
0+ = in the bulk liquid (i.e., g((f> s o) = g(0) = and g(<p) > otherwise). When boundary 
conditions are applied such that — » <f>± for z — > ±oo, a planar interface orthogonal to z is 
forced to appear in the system. The corresponding z-dependent OP profile is the extremal 
point 4>o(z) of (DO) that satisfies the boundary conditions: 

aj>Q - «0o" = ^(0o; T = T m ) , with </> (-oo) = (p s0 and (p (+oo) = . (2) 

Evidently, G[(f>o] represents the free-energy cost of the solid-liquid interface. 

Away from coexistence, the absolute minimum of g(<fi) falls at = <p s > for AT = 
T — T m < 0. This can be described by 

g{4>) = c 2 (j) 2 + c 3 3 + c 4 4 + . . . (3) 

if we take c 2 = C2o+c' 20 AT (c 2 o, c' 20 > 0), all other c n coefficients being constant. Not far from 



T m , the OP profile of a spherical solid cluster of radius R 3> a/2c/c2o is well described by 
0o (f — R), provided that the center of (po( z ) is chosen at z — 0. Based on these assumptions, 
the free energy of cluster formation becomes [2]: 



f+oa r+00 

AG(R) = 4tt / drr 2 [c0 o 2 (r - R) + 2/t0 o ,2 (r - R)] - 4tt4 |AT| / drr 2 2 (r - R) . 
Jo Jo 

(4) 
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A straightforward calculation then shows that: 

/ + OO 1 
dr r 2 2 (r - R) = -<p 2 R 3 + (l' 2 - J 2 ) + 2{I[ + h)R + (l' - I )R 2 ; 

n+OO 

/ dr r 2 0' o 2 (r - R) = J 2 + 2J X R + J R 2 ; 
Jo 

/■+00 

/ dr r 2 o /2 (r - i?) = K 2 + 27^7? + K R 2 , (5) 
where, for n — 0,1, 2: 

" + 00 /"+00 



n+oo r+oo 

In = / d^(0 2 o -0 2 (-z)) and /;= / dz 
Jo Jo 



z n <J> 2 (z) ; 



-+oo 

J„ = / dzz n (f)' 2 (z); 

oo 

+ 0O 



K n = / d^> /2 (z). (6) 

oo 



Substituting Eqs. (jSJ) into (Hj), we obtain the final expression for AG(R) (Eq. (2) of MT, with 
(7 = 0): 

AG(i?) = 4^ ^1 - ^+ -^j - ^vr^lA/xl , (7) 

where 

a = cJ + - 4 (/£ - J )|AT| ; 

aS = -cJ x - 2kK x + + Ji)|AT| ; 

ae = cJ 2 + 2kK 2 - c' 20 (I' 2 - h)\AT\ ; 

p s Ap = c'^AT . (8) 
In particular, the solid-liquid interface tension a and the "Tolman length" 5 are given by: 

/+oo 
dz [c<f)' 2 {z) + 2k^ 2 {z)] ; 
-oo 

1 f + ™dzz[c<P' 2 (z) + 2K<P" 2 (z)] 

Note that, for k = 0, the formula for 5 reduces to that reported in |3j. A nonzero 5 
corresponds to a <po(z) that is not symmetric around 0, namely to an interface between 
phases of a different nature. Summing up, Eq. (J7J) describes the corrections to the classical 
nucleation theory (CNT) which arise by replacing the assumption of a sharp solid-liquid 
interface with a more realistic finite width. 
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II. CALCULATION OF THE SURFACE TENSION AND THE TOLMAN 
LENGTH 



Given the form of g, one can compute the explicit values of a and 5 in Eqs. ([H]), and of 
e = (cJ 2 + 2kK 2 )/<j (i.e., the value of ?at T m , see Eq. (jSJ)), once the exact <j>o(z) is known. 
In turn, (f>o(z) follows from solving the boundary value problem (J2J), which can be simplified, 
after an integration by parts, to 

= + f 0(f - <7(0o) , with o (-oo) = sO and o (+oo) = . (10) 
A special case of g function is used in the MT, where at T m we take 

g(4>) = c 2 o^ 2 (l - j-X (l + t£-) = g (4>) (l + r£-) . (U) 



For this g, the differential equation (110]) is still too difficult to solve in closed form for generic 
k, even for r = 0. Hence, we decide to work perturbatively in k and r. 

At zeroth order, i.e., k = r = 0, corresponding to 4 theory, the solution to (ITU1) is 



0o^ = ^|l-tanh^^jj (12) 

with £ = a/2c/c2o and arbitrary C. We fix C by requiring that the interface is centered at 
z — 0, i.e., by imposing 

»+oo 

dzz0o(z) = O (13) 

(hence C = 0). Next, we take non-zero k and r, assumed to be of the same order of 
magnitude, and search for a first-order solution to ffTU]) in the form 

</>o{z) =fo{z)+Til) l (z) + ^xx{z)- (14) 

Upon inserting this function into Eq. (|T0|) . we obtain two independent equations for ipi{z) 
and XiW; namely 

c m-9m^=^f^ as) 

and 

<x'i ~ 9o(h)Xi = of (M - ^'o' 2 ) • (16) 
The solutions to Eqs. fll5p and (TT^j) such that each term in Eq. (1141) separately meets a 
requirement analog to ([TBI are: 

^■-iJUt'-^i-^i) (17) 
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and 

Upon plugging the by now specified (f>o(z) in the integrals defining a, 5, and e, we eventually 
obtain the final expressions quoted in Eq. (4) of MT. 



III. DERIVATION OF THE INTERFACE HAMILTONIAN 

We synthetically show how the interface Hamiltonian, Eq. (5) of MT, can be obtained 
from the Landau free energy ([1]). For this derivation, we build on Refs. % |5). More details 
will be given in [2|. 

In the same spirit of a statistical field theory for vesicles, we wish to assign a free energy 
cost to each particular realization of the cluster interface, here assumed to be sharp and akin 
to a mathematical surface E. The goal here is to describe the effect of shape fluctuations 
- though, in practice, in order to make analytical progress, we shall be forced in the next 
paragraph to describe just small deviations around the spherical shape. 

Let R(u, v) be the parametrization ("coordinate patch") of a tiny piece of E. For points 
r close to this small portion of E, we switch from 3D cartesian coordinates, r = (x, y, z), to 
new coordinates q a = (u,v,() (tangential and normal to E): 

r = R(u, v) + (n(u, v) , (19) 

where 

s(M= ijHnri (20) 

| -tin A ttv | 

is the unit normal to E. For a patch that deviates only slightly from planarity, we may 
adopt a free energy of Q[4> (((x,y, z))], thus arriving, by Eq.(JT]), to the surface Hamiltonian 

U s = J dnd^dC|J|{^(V0o(C)) 2 + ^(V 2 0o(C)) 2 + ^o(C))} (21) 

with J = \r u ■ (r v A r^)| = n ■ (r u A r v ). In order to compute the explicit form of J as well 
as those of the gradient and Laplacian of a function of £ only, it is convenient to take a 
patch parametrization in terms of orthonormal, arc-length coordinates, i.e., one such that 
R„ • R\, = all over the patch and |R U | = |R„| = 1. Although this can rigorously be done 
only for surfaces with zero Gaussian curvature (K = 0) [6], we can reasonably expect to 
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make small errors of order K for quasi-planar interfaces. It then follows that 

dr ^ 

dr 



n : 



— = (l - C4 1} )R« - (t 9 R v 



dr 

dv 



-Cr 3 R u + (l-C4 2) )R, 



(22) 



where Kn ^ and Kn' are the normal curvatures of the u- and w-lines respectively, and r g 



-(2) 



.(1) 



Tg is the geodetic torsion. From Eqs. (122 p . one can derive the metric tensor g a s of 



the transformation (I19p and the Jacobian: 



9a8 = 



dr dr 



dq a dq p 



( (l - (k^) + Cr 2 g ~2(r g + (\ + K l 2) ) \ 
-2(r 9 + (\ («« + ^ 2) ) (l - C^ 2) ) 2 + C% 2 



(23) 



and 

J=(l-C4 1) )( 1 -C*? ) )-C 2 ^ = ^ ) (24) 
g being the determinant of ( J23l) . Considering that covariant and contravariant components 
of a vector v = X]i=i ( v ' are built by projecting v on the bases Vq a and dr/dq a , 

respectively, we can calculate the gradient of a scalar field <fi and the divergence of a vector 
field A in local coordinates as follows: 



V0 



-9 



ad 



Or 



dq a dq/3 

gai3 Dem g inverse of f )23|) . In particular 



and V • A 



J__d_ 

y/gdq c 



(25) 



V0(C) = 0'(C)fi and V 2 0(C)=0 ,, (C) + 0'(C)V-n, 



(26) 



where 

V • n = i (-«?> - 4 2) " 2Cr|) . 
Finally, the mean and Gaussian curvatures of the patch are given by 



H = -V n 

2 



C=o 



2 i K n^ + K n 



and 



_ . dn <9n 
AT= n- ( — A — 

du dv 



C=o 



K (1) K (2) _ 2 
"'n ""n 'g ' 



(27) 



(28) 



(29) 



15 



We are now in a position to simplify Eq. (12 ip . Upon using Eq. (J2J) to eliminate g(<fro) in 
favor of (c/2)#£ - k [(^q)' ~ (3/2)#) /2 ], and inserting Eqs. ([21]), ([26]), (J2HD, and ([29]), we 
eventually get |2| the classic Canham-Helfrich Hamiltonian for fluid membranes, 

H a \E] = J dS (a + bH + cH 2 + dK) , (30) 

with the following explicit expressions for the coefficients: 

dC M 2 (0 + 2k$ 2 (0] ; 

oo 

6 = 2/ dCC[c0' o 2 (C) + 2^ 2 (C)] ; 

>/ — oo 

/+oo 
dC0o 2 (C); 
-oo 

/■+OO 

<* = / dC{C 2 [c0' o 2 (C) + 2^ 2 (C)]-^' o 2 (C)}. (31) 



oo 

/• + 00 



A few remarks are now in order: a) in deriving Eq. (1301) . all subleading corrections to the H 2 
and K terms were ignored, b) H and K are reparametrization invariants, hence no ambiguity 
arises from the arbitrariness of the parametrization used, c) The above derivation actually 
applies for just one S patch. However, upon viewing E as the union of many disjoint 
patches, the Hamiltonian (1301) holds for S as well, d) As anticipated, the coefficient d of 
the K term in ( 1301) is possibly different from the quoted one since a parametrization in 
terms of orthonormal coordinates does not generally exist. However, as far as we only allow 
for clusters with the topology of a sphere, J s dS K takes the constant value of 4ir by the 
Gauss-Bonnet theorem and the K term in % s can be dropped. Comparing the definition of 
a and 6 in Eqs. ([3T]) with Eq. ([9]), we can rewrite Eq. ( 130]) in the form 

H S \E] = [ dS (a - 2aSH + 2\H 2 ) , (32) 

where A = c/2. e) The term linear in H is related to the spontaneous curvature of S, 
H Q = — 6/(2c), which is proportional to the Tolman length 5. A nonzero value of H entails 
a difference in energy between inward and outward interface protrusions, with the effect of 
producing a nonzero Tolman length 5. The realization that in systems where, as in the Ising 
model, the symmetry is perfect between the two phases then 8 = 0, has long been known [3J. 
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IV. FIELD THEORY FOR THE INTERFACE HAMILTONIAN 



We consider a single solid cluster in a supercooled-liquid host. We model the cluster 
surface £ as being subject to random fluctuations distributed according to exp{— /3"H S [E]} 
with j3 = l/(k B T) and "H S [S] given as in Eq. (5) of MT. We stress that this form of interface 
Hamiltonian is only valid slightly below T m , where the curvature of the cluster surface is 
small. The cost in free energy of a cluster of volume V is then taken to be 

AG(V) = -p a \AfjL\V + F a (y), (33) 

where the surface free energy F s = —(1/(5) \nZ s with 

Z s (V)=a 3 J £>E5(V[£]- V)exp{-FH,\Z]} ■ (34) 

— 1 /3 

In the above equation, a = p s is a characteristic system length, V[S] is the volume enclosed 
by S, and £>£ is a suitable integral measure (see below). Neglecting overhangs and liquid 
inclusions in the cluster, let r = R(0, 4>) be the equation of £ in spherical coordinates. To 
proceed further, we assume only small deviations from a sphere, i.e., R(9, <fi) = Ro[l+e(9, (/>)], 
with e(6*,0) <C 1. Then, we expand e(#,0) in real spherical harmonics, 

oo I 

< 9 ^) = J2Y1 Xl,mYl, m {QA), (35) 
1=1 m=—l 

and we agree to ignore, from now on, all terms beyond second-order in the coefficients x^ m . 
With these specifications, the enclosed volume and area of £ are approximated as 

V[E] = ^kRI + R 3 J2 x lm = /({*}) (36) 

l>0,m 

and 

A\S[ = ^Rl + ^ ]T (/ 2 + Z + 2) x lm = 4 ^o <?({*}) , (37) 

l>0,m 

respectively, f({x}) and g({x}) being close-to-1 factors. In order to evaluate the mean 
curvature H, we move from 

V • n = (l + h 2 e(9, 0) - le(0, <P)L 2 e(9, 0)) , (38) 

where 

1 9/ J\ Id 2 



sin ^ -zns- ( 39 ) 



sin 9 09 \ 89 J sin 2 9 dcj) 2 
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Eventually, we obtain j^j: 

/ dS(a- 2a5H + 2XH 2 ) = AnaR 2 + ^ V(I 2 + I + 2) 
/v 2 ^ 



r 2 



l>0,m 

-braSRo - a5R ^ 1(1 + \)x\ m + 8n\ + ^ + *)(' " *)(' + 2 K™ • (40) 

l>0,m l>l,m 



Finally, we specify the integral measure in ( 134")) : 



00 *>0,m VS 7 70 



where S = (36vr) 1 / 3 1/ 2 / 3 is the area of the spherical surface of volume V and s = Ana . 
Equation ( I4ip follows from requiring that the present theory be the continuum limit of the 
field theory for a solid-on-solid (SOS) model whose heights are defined on nodes uniformly 
placed over a sphere of radius a/ S/ (Air) {2J]. 

We can now compute the partition function (1341) . We calculate first the integral on R 
by rearranging the delta function in Z s as 

* (girflS /({«}) -r)= v ww , A( , }r ; - <«> 

After doing the trivial integral over i? 0) we remain with a factor /({x})^ 1 ^ 3 which, within a 
quadratic theory, can be treated as follows: 

\ -1/3 

V .T? 1 

An 



f({x})^ 3 = ( 1 + A £ ^ J) exp j--L £ zjl . (43) 

l>0,m / l>0,m L Z>0,m J 



Z>0,m 

In the end, we arrive at a Gaussian integral which is readily computed: 

-2/3 f /,t/\ 1/3' 



Z s = (36tt)- 1/3 (J^J exp I Pp a \Afi\V - f3aS - 8n/3X + 8nf3a5 (J^J 

n 



2nS x 3 



z>i 



4vr/3(x5W^(Z 2 + / - 2) 
V An 



f3aS 2 



1 + ^-(l +l~2) + 2n/3X 1(1 + l)(l- l)(l + 2) 

-(i+l/2) 

(44) 



Without a proper ultraviolet cutoff Z max , the Z sum in F s = — IcbT In Z s does not converge. 
This is a typical occurrence for field theories on the continuum, which do not take into 
account the granularity of matter at the most fundamental level. We fix Z max by requiring 
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that the total number of (I, m) modes be equal to the average number of SOS heights/atoms 
on the cluster surface. It thus follows: 

S 



1 . 



(45) 



With this cutoff, the surface free energy becomes F s = 7(5)5, with a surface tension 7(5') 
dressed by thermal fluctuations: 

7(5) 



25 



( 2l + !) ln [ A + B (l 2 + 1 - 2 ) + C(/ 2 + 1 ~ 2 ) 2 ] 



1=2 
1/2 



SJ S \a 2 J S 

The quantities A, B, and C in Eq. (j4"6"l) are given by 

2C 



+ 



5 



A = ^ B 

5 2 ' 5 2 ' 5x/5 ' 5 



-Do . B Q _ Co 



where 



An 



-Dfl — _ „ , — — , Uq 



0a6s* 



4tt 2 ' ~ u 8tt 2 ' ~ u 2tt ' ~ u 2vr^' 



(46) 



(47) 



(48) 



By the Euler-Mac Laurin formula, the residual sum in Eq. ( 146]) can be evaluated explicitly. 
After a tedious and rather lengthy derivation, we obtain (for A 7^ 0): 



7(5) = a + 



k B T 
2a 2 



ln 



Bp 

a 2 e 2 



1 + ^lh.l I 



Co 



+ 
+ 



-2^+M^lnll 

4C O 0F 



Co 
5na 2 



1/2 



47TV- 7 ^ ln(5/a 2 ) 
"5 J "6^ T ^^" 



8^A-31n^-Hln^ + 3-?ln2-^ + 121 



D a 



D 2 



6 



96 46080 



4C 4B Co 6""V a2 fl4 / 8Co(5 a 2 + C ) 2 
x I -iB C D a 3 - ISfioCo'a 2 - 2C^D a - yC 3 - y5 2 Coa 4 



-2BlD a 5 + 2B Dla A + 2C Dla 2 )] 



S 



(49) 



up to terms o(5 _1 ). We wrote a computer code to evaluate the sum in (1461) numerically, and 
so checked that every single term in Eq. (14*91) is correct. In the notation of MT, the interface 
free energy 7(5) has the form: 

25 



7 = 0-1 



R + R 2 



7 JnjR/a) 



(50) 
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with R = [3V/(47r)] 1//3 and T-dependent expressions for a, 5, and ? that can be read in 
Eq. (149 p . At coexistence, the values of a and 5 for a quasi-spherical cluster differ from those 
in Landau theory (i.e., a and 8) for fluctuation corrections which are the effect of thermally 
excited capillary waves on the cluster surface. 
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